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WR20a: a massive cornerstone binary system comprising two 

extreme early-type stars* 
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Abstract. We analyse spectroscopic observations of WR20a revealing that this star is a massive early-type binary system 
with a most probable orbital period of ~ 3.675 days. Our spectra indicate that both components are most likely of WN6ha or 
03If*/WN6ha spectral type. The orbital solution for a period of 3.675 days yields extremely large minimum masses of 70.7±4.0 
and 68.8 ± 3.8 M G for the two stars. These properties make WR 20a a cornerstone system for the study of massive star evolution. 
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1. Introduction 

Although it is generally accepted that Wolf-Rayet (WR) stars 
are evolved descendants of massive O-stars, many open ques- 
tions remain about the evolution of these objects. Among the 
most fundamental issues is the importance of mass transfer in 
the evolution of massive O + O binaries (e.g. Vanbeveren et 
al. I1998I . It is therefore of prime interest to achieve as com- 
plete a census of WR binaries as possible and to investigate the 
fundamental properties of individual massive binary systems. 

Among the 227 objects listed in the Vllth Catalogue of 
Galactic Wolf-Rayet Stars (van der Hucht 2001 1, there are only 
52 known spectroscopic binaries and only 32 of them have re- 
liable orbital solutions. In addition to these confirmed bina- 
ries, van der Hucht (2001 1 provides a list of 30 rather faint 
and poorly studied objects that are potential WR + OB bina- 
ries. Most of these binary candidates were classified as such 
because of their comparatively weak emission lines that could 
be the result of the dilution by the continuum light of an OB- 
type companion. In order to search for a composite spectrum - 
and hence confirm multiplicity - we performed a snapshot spec- 
troscopic survey of a sample of WR binary candidates. One of 



the most interesting objects that emerged from our survey was 
WR20a. 

WR20a (= THA 35-H-036 = SMSP2) was catalogued by 
The (1966) in an Ha emission line survey, and first discov- 
ered to be a WR star by Shara et al. ( 1991 ) in the course of 
a deep narrow and broadband photometric survey for galactic 
Wolf-Rayet stars in the southern Milky Way. The latter authors 
also obtained a low -resolution (2 pixel resolution of 7 A) spec- 
trum of the star and proposed a WN7 spectral type, although 
they noted the lack of a He n A 541 1 emission line. Shara et al. 
dl991t suggested that WR 20a had a rather high H/He abun- 
dance ratio. Later, the object was re-classified as WN7:h/WC 
by Shara et al. J1999I . The star is believed to be a member of 
the young open cluster Westerlund 2 in the core of the H n re- 
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gion RCW49 (Moffat et al. HWTt . 

In this letter, we first present our blue/violet observations of 
WR20a (Sect.|2j. The spectral features as well as the nature of 
WR 20a are revisited in Sect. 13.11 whilst Sect. l3.2l deals with an 
analysis of the radial velocities. Finally, we briefly discuss the 
implications of our results in Sect.0] 

2. Observations 

Two snapshot spectra of WR 20a were obtained in March 2002 
with the EMMI instrument mounted on ESO's 3.5 m New 
Technology Telescope (NTT) at La Silla. The EMMI instru- 
ment was used in the BLMD medium dispersion spectroscopic 
mode with grating # 3 (1200 lines mm -1 ) providing a wave- 
length coverage of 460 A per set-up. The central wavelength 
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Fig.l. Spectrum of WR20a as observed on HJD 2453 034.715 (violet setting, top panel) and 2453 034.741 (blue setting, lower 
panel). Note the double N iv A 4058 emission line as well as the double N v A 4604 and N v A 4620 absorptions. The stellar lines 
are labelled above the spectrum, whereas the interstellar lines and diffuse interstellar bands (DIBs) are indicated below. 



was set to 4080 A (violet setting) or 4530 A (blue setting). The 
exposure times were 30min for each setting and the corre- 
sponding S/N was about 100. The spectral resolution as de- 
termined from the FWHM of the lines of the Th-Ar comparison 
spectra was about 1 .25 A. The data were reduced in the stan- 
dard way using the long context of the midas package. The 
violet spectrum taken on HJD 2 452 355.60 revealed double ab- 
sorption lines of hydrogen and helium as well as a moderately 
strong double-peaked Niv A 4058 emission. On the other hand, 
the blue spectrum obtained on HJD 2 452 354.56 showed single 
absorption lines of hydrogen, He n and N v. These properties 
suggested that WR 20a could be a short period binary consist- 
ing of two very early extreme Of-type stars (Driesens 2003 1. 
We thus organized a follow-up observing campaign to confirm 
the binarity of WR20a. The star was monitored over 12 con- 
secutive nights in January-February 2004 with the same instru- 
mentation as in March 2002. The weather conditions were good 
except for the last two nights where the seeing and the trans- 
parency were poor. The data were reduced in the same way as 
in 2002 and the complete journal of the useable observations is 
given in Tabled 

3. Results 

3.1. The spectrum of WR20a 

Our data reveal that WR20a is indeed a short-period binary 
system (see below). All prominent spectral features appear 
double with nearly identical intensities on most of our spectra 



(see e.g. Fig.Q indicating that the components of WR 20a must 
have identical spectral types and a visible brightness ratio of 
~ 1 . The spectrum of WR 20a displays many absorption lines 
of H i, He n and N v unlike what we would expect if WR 20a 
were indeed a 'classical' WN7 star (see however below). Only 
a few (rather narrow) emission lines are seen in the blue/violet 
spectrum: N iv A 4058, He n A 4686 as well as a broad complex 
probably due to Nm AA 4634-41, Cm AA 4647-50 and possi- 
bly C iv A 4658. The Balmer absorption lines also display some 
broad underlying emission wings on several spectra. The emis- 
sions seen around the H6 absorption in Fig.^could actually be 
due to Si iv AA 4089, 4116, but the spectra obtained on other 
nights do not allow us to unambiguously confirm this identifi- 
cation. In addition to the strong diffuse interstellar bands, we 
note the presence of the interstellar He i A 3889 absorption line. 
This latter feature was also observed in the spectra of several 
early-type stars of the Carina Nebula (Walborn & Hesser [T975l . 
This line is not usually seen in interstellar sightlines (cf. £ Oph, 
Shulman et al. 1 19741 1 and may be formed locally as in the case 
of the Carina stars. 

No stellar Hei lines are clearly seen in the spectrum. 
Together with the presence of rather strong Nv absorption 
lines, this would suggest a very early O-type classification 
for the components of WR20a. The presence of Hen A 4686 
in strong emission further suggests a supergiant luminosity 
class. Comparing our spectra to the classification scheme of 
Walborn et al. (2002 see also Walborn 2001 1, we would pro- 
pose an 03 If* spectral type for both stars. However, distin- 
guishing an extreme Of star from a weak-lined WNL star 
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is a very tricky issue. In this context, Crowther & Dessart 
l fT998l proposed EW(Hen A 4686) = 12 A as a boundary be- 
tween 03 If*/WN and WN5-6 stars. For WR 20a, we determine 
EW(Hen ,14686) = 13.1 A with a 1-cr dispersion of 1.3 A, 
right at the borderline proposed by Crowther & Dessart. Since 
the He u A 4686 emission could contain a contribution from a 
wind interaction process, we consider this an upper limit on the 
actual strength of the line. Comparing the spectrum of WR 20a 
obtained near conjunction with the spectra of several extreme 
Of and weak-lined WN stars, we find that the components of 
WR20a are indeed intermediate between Mk42 (03If/WN) 
and HD 97950C (= WR47c, WN6 ha). Fi nally, the spectrum of 
WR 20a presented by Shara et al. J199U indicates that the U/3 
line is in weak emission. The latter feature favours a WN6ha 
spectral type. We shall thus adopt a WN6ha classification for 
both components of WR 20a, though we cannot fully exclude 
an 03 If /WN6ha type. The Shara et al. JTWlT l spectrum of 
WR 20a suggested the presence of a C iv A 5808 emission with 
a normalized intensity stronger relative to Hei A 5876 than in 
normal WN stars, thus leading to the WN/WC classification 
proposed by Shara et al. ( 1999 ). Note that our blue/violet spec- 
tra do not provide support for such a hybrid type. 

3.2. Orbital solution 

We have measured the radial velocities (RVs) of the compo- 
nents of WR20a by fitting Gaussians. The typical uncertainties 
on individual RV measurements are of order 10 - 15 km s _1 for 
those phases when the components are resolved 1 . 

A major issue in determining the orbital period arises from 
the fact that the spectral types of both stars are identical. In 
fact, the spectra do not allow us to distinguish the components 
a priori. Therefore, we performed a period search on the abso- 
lute values of the radial velocity differences \RV\ — i?V"2l- In this 
way, we expect to find the highest peak in the periodogram at 
vi = 2/f 01 b. Using the Fourier technique of Heck et al. ( 119851 
see also Gosset et al. 2001 1, we find that the most likely value 
of vi is 0. 544+0. 009 day -1 . However, there is a strong ambigu- 
ity with the 1 — vi alias at 0.453 + 0.009 day -1 which produces 
only a slightly smaller (formal) peak in the periodogram (see 
Fig.fJ}. Our preferred value for the orbital period is therefore 
3.675 + 0.030 days, but we cannot discard the alternative value 
of 4.419 + 0.044 days. Alternatively, if the (slightly smaller) 
peak at the 1 + v\ alias (1.547day~') were to be selected, the 
orbital period would be ~ 1.29 days, i.e. even shorter than the 
1 .64 day period of CQ Cep, the shortest period WN binary sys- 
tem known to date (see e.g. van der Hucht 2001 1. 

As a next step, we determined the apparent systemic veloc- 
ity of each spectral line from an orbital solution obtained for the 
RVs of the specific line. In this way, we found apparent sys- 
temic velocities of -75, -18, -77, -42, -51 and -26kms~ 1 
for the Hy, Hen ,14542, Nv ,14604, Nv ,14620, Niv ,14058 

1 Note that we also determined the heliocentric radial velocities of 
the sharp interstellar features: we obtain RVs of -16.6 + 15.0, -7.8 ± 
10.5 and -10.2 ± 14.0 kms" 1 for the Hei ,13889, Can ,13934 and 
Can A 3968 lines. The lcr dispersions listed here provide some rough 
indication of the accuracy of our wavelength calibration. 
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Fig. 2. Top panel: power spectrum of our \KV\ - 7? V2I time se- 
ries of WR 20a. The bottom panel provides the corresponding 
power spectral window. 



and H<S lines respectively. Finally, we averaged the RVs of the 
different lines listed above after subtracting their systemic ve- 
locities. The results are listed in Tabled 

Due to the ambiguity about the right alias in the peri- 
odogram, we caution that the orbital phases as well as the pri- 
mary/secondary identifications in Tabled are only valid if the 
orbital period is indeed 3.675 days. Adopting this period, we 
have derived the orbital solution of WR 20a using the technique 
described by Sana et al. (2003 1. Since the RVs show no indica- 
tion of a significant orbital eccentricity, we assumed a circular 
orbit. The results are provided in Tableland the RV curve is il- 
lustrated in Fig. [5] We emphasize that orbital solutions derived 
from the RVs of individual spectral lines overlap within the er- 
rors with the results in Table|2] 

The most important results are the large minimum masses 
of 68.8 + 3.8 and 70.7 + 4.0 M Q derived for the secondary and 
primary component respectively. It should be emphasized that 
these values strongly depend on the adopted orbital period. In 
fact, the alternative orbital period (4.419 days) would lead to 
even larger minimum masses of 80.8 ±5.1 and 84.0 + 5.3 M 
respectively. If, on the contrary, the 1 + v\ alias corresponding 
to an orbital period of 1 .293 days were selected, the minimum 
masses would be much lower (24.9 +1.4 and 25.9 ± 1.4 M ). 
We note however, that in the latter case, the orbital solution is of 
significantly lower quality than the one obtained for an orbital 
period of 3.675 days. 

4. Discussion and future perspectives 

The WN6ha or 03 If*/WN6ha spectral type of the components 
of WR20a is significantly earlier than those of the earliest 
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Table 1. Journal of our spectroscopic observations of WR 20a. 
The first and second columns provide the heliocentric Julian 
date (in the format HJD-2450000) and the orbital phase as- 
suming a period of 3.675 days and adopting phase <p - 0.0 as 
given in Tabled The third column indicates whether the blue 
(B) or violet (V) spectral range was observed. The recentred 
and averaged primary and secondary radial velocities are given 
in columns 4 and 5 respectively (see text). Finally, the last col- 
umn yields the weight that we assigned to the various RV data 
in the orbital solution. 



Date 


<P 


Set-up 


RV] 
(kms- 1 ) 


RV 2 

(kms" 1 ) 


weight 


2354.560 


0.542 


B 


-11.7 


-11.7 


0.1 


2355.600 


0.825 


V 


290.7 


-319.3 


0.7 


3032.716 


0.074 


B 


-131.9 


182.9 


0.2 


3032.742 


0.081 


V 


9.1 


9.1 


0.1 


3033.715 


0.346 


V 


-288.4 


277.4 


0.7 


3033.739 


0.353 


B 


-285.1 


292.6 


1.0 


3034.715 


0.618 


V 


235.8 


-219.1 


0.7 


3034.741 


0.625 


B 


287.6 


-265.5 


1.0 


3035.719 


0.892 


V 


269.3 


-222.7 


0.7 


3035.746 


0.899 


B 


245.2 


-219.4 


1.0 


3036.719 


0.164 


V 


-259.7 


287.0 


0.7 


3036.743 


0.170 


B 


-302.6 


309.5 


1.0 


3037.738 


0.441 


V 


-78.3 


-78.3 


0.1 


3037.762 


0.447 


B 


-74.0 


-74.0 


0.1 


3038.723 


0.709 


V 


363.5 


-354.0 


0.7 


3038.746 


0.715 


B 


374.6 


-375.9 


1.0 


3039.718 


0.980 


V 


12.4 


12.4 


0.1 


3039.752 


0.989 


B 


-28.0 


-28.0 


0.1 


3040.717 


0.252 


V 


-309.3 


368.4 


0.7 


3040.741 


0.258 


B 


-336.1 


394.0 


1.0 


3041.722 


0.525 


V 


4.3 


4.3 


0.1 


3041.747 


0.532 


B 


22.9 


22.9 


0.1 


3042.714 


0.795 


B 


375.4 


-374.4 


0.2 



Table 2. Orbital solution for WR 20a assuming a circular orbit 
and an orbital period of 3.675 days. To refers to the time of 
conjunction with the primary being behind, y, K and a sin i 
denote respectively the systemic velocity, the amplitude of the 
radial velocity curve and the projected separation between the 
centre of the star and the centre of mass of the binary system. 
Rrl stands for the radius of a sphere with a volume equal to 
that of the Roche lobe computed according to the formula of 
Eggleton ( 1983). All error bars indicate 1-cr uncertainties. 



Primary Secondary 



T (HJD-2450000) 


3043.480 + 0.014 




y (kms -1 ) 


18.8 + 5.6 -3.5 ± 


5.7 


if (kms" 1 ) 


353.1 +8.6 362.6 ± 


8.8 


a sin (Rq) 


25.6 ±0.2 26.3 + 


0.6 


q = m l /m 2 


1.03 + 0.03 




m sin 3 (M Q ) 


70.7 ±4.0 68.8 ± 


3.8 


Rrl sin (' (R Q ) 


19.8 ±0.2 19.5 ± 


0.2 
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members (06:V - 07: V) of Westerlund 2 analysed by Moffat 
et al. fTWTY 

Several authors have estimated the distance of 
Westerlund 2 : 7.9! j;gkpc (Moffat et al. [T99TV . 5.7 ± 0.3 kpc 
(Piatti et al. fT998t : 6.4 + 0.4 kpc (Carraro & Munari I20041 
Adopting a visual brightness ratio of 1 together with the 
photometry of Moffat et al. ( IT99T1 V = 13.58, B -V = 1.51) 
the various distance estimates yield absolute magnitudes My 
in the range -5.0 to -6.1. This is fainter than the average of 
three WN6ha stars in NGC3603, namely -6.9 to -7.7 mag 2 , 
such that only the largest published distance of Westerlund 2 
would be in marginal agreement with the photometry of 
WR20a, allowing for a cosmic scatter of 0.5 - 0.7 mag 
on My. This problem would be less severe if the absolute 
magnitudes of the components of WR 20a were more typical 
of 031 stars (My = -5.9 according to Crowther & Dessart 
1998 1. Alternatively, WR20a could actually be unrelated to 
Westerlund 2 (though the rarity of early-type stars and the 

2 Moffat et al. (E>94 1 derived an average My of —7.7 for the three 
H-rich WN stars in NGC3603, whilst Crowther & Dessart (1998 1 
re-evaluated their distance and reddening and obtained an average of 
-6.9 mag. 



Fig. 3. Radial velocity curve of WR 20a for an orbital period 
of 3.675 days. The filled (resp. open) symbols stand for the 
primary (resp. secondary) RVs. The sizes of the symbols re- 
flect the relative weight assigned to the various points in the 
orbital solution (see TableQ. The solid line and the dashed line 
yield our best fit orbital solution assuming a circular orbit (see 
Table©. 



angular proximity of WR 20a to the cluster strongly suggest 
membership) or the photometry of Moffat et al. could be af- 
fected by photometric eclipses. Concerning the last possibility, 
we note that given the extremely large minimum masses of 
the components in WR20a, it appears very likely that the 
system could display eclipses. A photometric monitoring of 
this binary is therefore of the utmost importance in order to 
help derive its absolute magnitude, to get rid of the aliasing 
problem of the RV data and to constrain the orbital inclination. 
We intend to organize such a campaign in the near future. 

An important issue in the context of the formation and the 
evolution of very massive stars is that of the mass of the most 
massive ones. Masses up to ~ 200 M have been inferred for 
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the most luminous 02 stars from a comparison with evolution- 
ary tracks (Walborn et al. 2002 1. However, to date, none of the 
02 stars is known to belong to a binary system that would allow 
to determine its mass in a less model dependent way. The most 
massive main-sequence stars in binary systems known so far 
were found in R 136-38 (Q3V + 06 V) and R 136-42 (03 V 
+ 03 V). Massey et al. (2002) determined masses of 57 and 
40 M respectively for the 03 V primary stars of these systems. 
Another extremely massive object is the WN7ha primary in the 
80-day period binary WR 22 (Rauw et al. 1996 Schweickhardt 
et al. 1999 1 with a minimum mass of 55 M Q . WR 20a thus prob- 
ably consists of two of the most massive stars with a direct 
mass determination known so far. Therefore, WR 20a is cer- 
tainly a cornerstone system for future investigations of massive 
star evolution. 

Finally, we note that the emission lines in the spectrum of 
WR 20a display strong phase-locked profile variability. This as- 
pect will be discussed along with an analysis of the spectrum 
with a model atmosphere code in a forthcoming paper. 
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